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Abstract: We present a method for the direct extraction of rational contributions to 
one-loop scattering amplitudes, missed by standard four-dimensional unitarity techniques. 
We use generalised unitarity in Z) = 4 — 2e dimensions to write the loop amplitudes in 
terms of products of massive tree amplitudes. We find that the rational terms in 4 — 2e 
dimensions can be determined from quadruple, triple and double cuts without the need for 
independent pentagon contributions using a massive integral basis. The additional mass- 
dependent integral coefficients may then be extracted from the large mass limit which can 
be performed analytically or numerically. We check the method by computing the rational 
parts of all gluon helicity amplitudes with up to six external legs. We also present a simple 
application to amplitudes with external massless fermions. 
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1. Introduction 



The main goal of the Large Hadron Colhder, due to start experiments later this year, is to 
explore the electroweak symmetry-breaking scale and search for new physical phenomena 
at the TeV scale. In order to effectively achieve this goal it will be necessary to have precise 
predictions of backgrounds within the Standard Model and particularly for the enormous 
number of QCD or QCD-associated events that could mask the effects of new particles. 

At the present time much progress has been made towards completing calculations 
for important cross sections [1]. For a selected number of processes, NNLO precision 
may be necessary but for the majority of events it is expected that NLO predictions will be 
sufficient. Because generic signals of new phenomena will be associated with the production 
and subsequent decay of heavy particles, the main source of backgrounds comes from 
multi-jet final states. One of the most important ingredients for such multiparticle NLO 
cross sections is the virtual matrix elements. The traditional Feynman approach to such 
calculations is extremely computationally intensive due to the rapid growth in the number 
of diagrams with the number of external legs. 

On-shell techniques, pioneered by Bern, Dixon and Kosower during the mid-nineties 
[2, 3], offer an elegant alternative to the traditional Feynman approach. When constructing 
loop amplitudes from on-shell objects one works only with the physical degrees of freedom 
which can substantially reduce the complexity of the calculation. Loop amplitudes are re- 
constructed by sewing together products of on-shell tree amplitudes to extract information 
about the branch cuts in each channel of the momentum invariants. This information can 
then be used to find the coefficients of the known scalar integral basis [4-6]. In supersym- 
metric theories performing the on-shell cut in four dimensions is sufficient to reconstruct 
the full amplitude, whereas in QCD we will be left with additional rational terms. Knowl- 
edge of universal factorisation and use of triple as well as double cuts [7] can be used to 
supplement this approach and reconstruct full amplitudes in some cases to all multiplicity. 

More recently on-shell techniques have seen renewed interest, via excursions into 
twistor space, through developments exploiting the use of complex momenta. On-shell 
recursive techniques at tree-level [8, 9] use basic complex analysis and universal factori- 
sation properties to write simple relations between on-shell amplitudes. These relations 
have been used successfully to derive compact analytic expressions for a wide variety of 
processes. Since the use of complex momenta ensures that three-point amplitudes are well 
defined on-shell, it is also possible to refine the generalised unitarity procedure and iso- 
late individual coefficients of the integral basis [10]. Using spinor integration techniques 
[11, 12] it has been possible to derive analytic expressions for the cut-constructible parts 
of all six-point gluon amplitudes [11, 13]. Such techniques are not restricted to massless 
theories and have also been developed for massive theories. They can also be applied to 
D-dimcnsional cuts [14-20], with recent applications to some examples of complete gluon 
amplitudes [21, 22]. 
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Given that the integral coefficients are rational functions of the external momenta it 
makes sense to look for a purely algebraic procedure which avoids any explicit integration. 
Cutting four propagators in four dimensions completely freezes the loop integral which 
allows integration to be replaced by algebraic operations [10]. After integrating a generic 
triple cut over the on-shell delta functions a one dimensional integral still remains and 
it is difficult to separate any new independent functions from the residues of previously 
computed, higher order poles. Ossola, Papadopoulos and Pittau showed that by a special 
parametrisation of the loop momentum one can find a systematic way to compute these 
higher order terms [23]. After a subtraction of such terms at the integrand level the 
remaining terms can be found by solving an algebraic system of linear equations. This 
method has been shown to be successful in the context of Feynman diagram calculations 
[24-26] and have been implemented in a public code CutTools [27]. 

A powerful analytical approach is the one by Forde who used simple complex analysis 
to go beyond the OPP approach and isolate the coefficients of the scalar integrals [28]. 
Introducing a complex parametrisation for the loop momentum, the coefficients are then 
completely determined through the limiting behaviour of the products of tree amplitudes. 
In particular this avoids the need to solve an algebraic system of equations and leads directly 
to compact analytic expressions for the coefficients. One can understand this procedure 
further by observing that the non-trivial integrations are simply contour integrals in the 
complex plane which, after application of Cauchy's theorem, give compact descriptions 
of triangle coefficients [29]. This method has been generalised to accommodate arbitrary 
internal masses [30] and, with modifications to remove higher order poles from the complex 
plane, shown to be an efficient numerical tool [31]. The technique has also shown to be a 
powerful tool for analysing non-trivial cancellations in gravitational theories [32] and for 
analytic computations of multi-photon amplitudes [33, 34]. 

Two alternative approaches to the computation of the rational terms have been followed 
in recent years. The first of these is the use of loop-level recursion relations exploiting the 
multi-particle factorisation properties [35-37]. This method has been used successfully to 
derive analytic expressions for many helicity amplitudes up to eight final state gluons [38], 
as well as all-multiplicity expressions for one-loop MHV amplitudes [39, 40]. The method 
applies equally well to amplitudes with massive external particles where amplitudes with 
gluons coupling to a Higgs boson have also been derived [41-43]. A good review of these 
new techniques together can be found in reference [44]. These recursion relations have 
been combined with the coefficient extraction of Forde to produce an automated C++ code, 
BlackHat [31]. Computation of amplitudes with up to eight external gluons have shown 
promising speed and accuracy at fixed precision. 

The second method is to use /^-dimensional cutting techniques which also completely 
determine the loop amplitude [45, 46]. Recently much progress has been made towards 
a numerical approaches to these techniques, based on the OPP approach. This has been 
shown to be much faster [47] than the current Feynman based techniques [48-52]. Giele, 
Kunszt and Melnikov [53] have used an approach using higher integer dimensions to provide 
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a purely numerical procedure for the evaluation of the D-dimensional coefficients. This has 
been implemented in a Fortran code, Rocket, and used to compute gluon amplitudes with 
up to 20 external legs [54]. Very recently this method has also been applied to amplitudes 
with massive fermions [55]. Ossola, Papadopoulos and Pittau have also proposed a tech- 
nique to calculate these rational parts [26] using a massive integral basis which motivates 
the construction we present here. 

The purpose of this paper is to extend the current D-dimensional approaches to com- 
putations of rational terms by making use the complex analysis techniques used by Forde 
for the cut-constructible terms [28]. This allows us to compute analytic expressions for 
full one-loop amplitudes within a single framework and sheds further light on the efficient 
computations of full one-loop amplitudes. Just as in Forde's analysis, this avoids the need 
to solve an algebraic system of equations present in the OPP approach and leads directly 
to compact analytic expressions. The main simplification arising from this analysis is the 
separation of pentagon contributions, which vanish in the four-dimensional limit and so 
are eliminated from the computation. This leaves the rational terms in terms of tree am- 
plitudes evaluated in the large- mass limit of box, triangle and bubble cuts. Alternatively 
this procedure can be understood as a contour integration for a complex mass parameter 
where the radius of the contour is taken to infinity. For the main part of this paper we 
concentrate on amplitudes with external gluons though the methods presented should also 
apply to more general external states by using the full D-dimensional tree amplitudes. This 
is demonstrated using a simple example involving massless external fermions. In section ^ 
we review the Z)-dimensional integral basis and the general form of the rational terms. In 
section ^ we describe how each of the components of the 4 — 2e-dimensional basis can be 
determined from the large mass behaviour of four-dimensional, massive, generalised cuts. 
Section ^ we present some analytic examples of gluon amplitudes with up to six external 
legs and outline a simple numerical implementation. We then present a simple application 
to a four-point massless fermion amplitude in section ^ before giving our conclusions. Some 
additional notes on mass dependence of the one-loop integrands and general forms of the 
boundary expansions are given in an appendix. 

2. Notation 

Throughout this paper we will be considering colour-ordered helicity amplitudes. We use 
the standard spinor-helicity formalism to describe all momenta and external wavefunctions 
[56-61]. Two component Weyl spinors are written as 



The indices of the two-component spinors are raised and lowered using two totally anti- 
symmetric tensors ea/3 and e^^ where ei2 = = 1- The scattering amplitudes are then 
written in terms of spinor products defined by, 






A°(p)A„((?) = {pq) 



Ad(p)A"(g) = [pq]. 
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Throughout the paper all momenta are written in matrix form via contraction with the 
Pauli a matrices: 

p" = {pla^lp], p.aaa = \p)a\p\a, (2.3) 

where we use a shorthand notation, 

p=\p)[p\. (2.4) 

For numerical evaluation of the spinor products we have used the standard approach as 
outlined in reference [62]. Modifications to accommodate complex and massive momenta 
have recently been implemented in a Mathematica package [63]. 



3. D-dimensional cuts and rational terms 



In this section we review the Z?-dimensional integral basis as considered in previous con- 
structions [14, 26, 53]. Here we focus on the connection between the 4 — 2e dimensional 
representation and that associated with an effective mass ^u^ [45, 46, 64]. 

We begin by writing a general 1-loop amplitude in terms of a D-dimensional n-point 
function, 

" J (47r)^/2 (/2 _ ml){{l - K^Y -ml)... ((/ + K^Y - mlY ^ ' ^ 

The numerator function N contains all information from external polarisation states and 
wavefunctions and tensor structures from the loop momenta. Since we are concerned with 
computing the four dimensional limit it is useful to decompose the loop momenta as, 

r = P + /f_2,], (3.2) 

where I contains the four-dimensional components and ^[-2e] contains the remaining D—A = 
— 2e dimensional components. Using D-dimensional Passarino-Veltman reduction tech- 
niques on ( p. ID allows us to reduce to a basis of scalar integral functions with rational, but 
/^-dimensional, coefficients [53], 

+ Y,CzMD)lE,K, + Y.^2MD)I?,K,+Ci{D)I?, (3.3) 

K-i K2 

where we define the sets of external momenta, Kr, as the set of all ordered partitions of 
the n external particles into r distinct groups (the ordering is defined by that of the full 
amplitude An^). We proceed by writing the amplitude in terms of an integral basis with D 
independent coefficients at the cost of expanding the basis of integral functions. Working in 
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the four dimensional helicity (FDH) scheme^ [66, 67] we keep all the external momenta and 
sources in four dimensions which means that the dependence on D can only arise through 
contracting the loop momentum with itself, 

f = l' + ll_^,^^l'-l?. (3.4) 

We can now interpret all dimensional dependence of the coefficients in terms of their de- 
pendence on /i^ . We also know that in any renormalisable gauge theory the maximum rank 
of an n-point tensor integral appearing in the amplitude is n, hence for the box functions 
we can have up to a maximum power of /i^ in the coefficient and up to ^u^ in the trian- 
gles and bubbles. The pentagon integral is only an independent function in D dimensions 
since we can find poles in the D — A dimensional sub-space. As a result, the coefficient of 
this function in D = 4 — 2e, or residue around the extra dimensional poles, can have no 
dependence on e. Therefore we arrive at a new basis, 

+ E cfkAK. w + E ctkAK. [A + E (^^SkAk. \A + E w 

Ki Ki Ki K3 

+ E ^S^3^3% [/^'] + E C2;kJ^,k, [1] + E C^^A, [f^'] + C^I[>. (3.5) 

K3 K2 K2 

The integrals over /j? can be performed by separating the integration into 4 and D — A 
dimensional parts, 

d^'h _ f d-^if,^) f d% 

(27r)^ J (2^)-2^y (2^)4- ^""-^^ 

It is fairly straightforward to write the four new integrals in terms of higher-dimensional 
scalar integrals using [45, 46], 



i^i^^^i = ^i^^'^'m ii(D-^+k). (3.7) 

k=0 

We also use the dimensional shift identity [68] to decompose the pentagon integrals^ 



Sij = l{m^ + m]-pl). (3.9) 



""^It is straightforward to convert gluon amplitudes into the 't Hooft Veltman scheme by subtracting a 
factor of ^Atioc from the result in the FDH scheme. Similar relations exist for amplitudes with external 
fermions [65]. 

^We refer the reader the [68] for full definitions of the quantities in eq. (3.8) although we point out that 
all integrals include a factor of ( — 1)"'+^ in their definition. The essential information is that the coefficients 
are just functions of the external momenta and the internal masses. 
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In the above is one of the five sets of four partitions obtained cychcally merging two 
adjacent partitions of a given pentagon configuration K^. After this is done the exphcit 
D-dependence of the amphtude is restored [53]: 

L>-4 



2 5;K5 



+ 2^ (~^4;Kji;Ki H ^ 2^ '-'A;kJ4;K4 ~^ 4 *^4;ii-4-'4;ii-4 + 2^ '-'^\Kih;K3 

K4 Ki K4. K3 

+ ^ E ^^3^3%^ + E C2;Kj,% + ^ E ^a.^S.^ + ^lA^' (3-10) 

where 

i=l j 

After taking the 4-dimensional hmit, D 4— 2e, we find that the integral basis reduces to a 
combination of box, triangle and bubble integrals but at the cost of introducing additional 
rational terms, 

4'^''-'^ = E^4;i.44^; + E^3;i.34^^ + E^2;X./2'^; + Cl/^2^ + i^„^ (3.13) 

K3 K2 

The final step is to identify how the higher dimensional integrals in ( 3.10| ) contribute to 
the rational terms. This is surprisingly simple since the scalar box and scalar pentagon 
integrals are finite in 6 — 2e dimensions and don't contribute to the rational part. The 
remaining three terms, written in terms of integrals over fi, are. 



4-2er 4i 

^ 6' 



T4-2er . 



ltV]'^'-l{s-3{ml + ml)). (3.14) 
The rational terms are thus given by [18, 26, 53], 

= E ^S.. - ^ E - I E - 3(-? + -^)) C^, (3.15) 

K4, K'i K2 



4. Extracting the integral coefficients using massive propagators 



To extract the integral coefficients using generalised unitarity we need to solve the con- 
straints which put the various propagators on-shell [28]. To generalise this to the D- 
dimensional case we also need to extract the ^ dependence of the coefficients as defined 
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in eq. ( ^.51) . Since the extra dimensions in the loop momentum can be interpreted as 
an effective mass term, it is possible to construct the full amplitude from tree amplitudes 
where the internal legs have a uniform mass: 



If = I 



72 







If 



A* 



(4.1) 



This method has been used successfully within the standard unitarity cut technique [45, 46] 
and in conjunction with spinor integration [18, 21, 22]. Solving the system of on-shell 
constraints can then be achieved in exactly the same way as the four-dimensional massive 
case [30]. For the current study we also consider ourselves to be restricted to cases with 
D-dimensional scalars with massless external fermions and gauge bosons. 



4.1 Box Coefficients 



h 



K4 



k 



K2 



Figure 1: A general quadruple cut with loop momentum flowing clockwise and all external mo- 
menta outgoing. 



In this section we will show that by extracting the coefficient of the Z)-dimensional box 
directly using Forde's formalism [28], we can ignore the pentagon coefficients entirely. We 
begin by choosing the four-momentum, /i, to be parametrised by, 



/- = aK\ + hK\ + c\Ki)[K\\ + d\K\)[Kl\. 
where the K\ ^ define a massless basis in terms of two of the external momenta: 



(4.2) 



t _ 714 (7l4i^4 - ^4-^1) 714 {IIaKi - SiK^) K K ^ fTv jr\2 

-^^4 - 2 ETT; ' ^1 = 2 ' 714 = -t^l • -^-4 ± V (^1 • -^^4)^ - '5l'i'4, 

7^4 - 6164 7f4 - 

(4.3) 

for Si = Kl 

The four on-shell constraints then fix the coefficients: 

h = aKl + bK\ + c\Kl)[K\\ + ^-^^^^^\k\)[kI\ 



C714 



C714 



\K\)[K\ 



(4.4) 
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where 

5'i(S'4 + 7i4) , 54(51+714) -ci ± Y^cf - 4coC2 , . 

'^ = — ErT^> = 2 ETT^' c± = , 4.5) 

7i4 - •5i54 7f4 - 5i54 2C2 

C2 = {k\\K2\k\1 (4.6) 
ci = a{Kl\K2\Kl] + 6(i^J|K2|i^J] - ^2 - 2Ki • K2, (4.7) 

co= fa6-^VE:J|i^2|i^l]. (4.8) 

V 714/ 

For the quadruple cut we find that both solutions for 714 are degenerate. 
Now consider the quadruple cut: 



=(4vr)^/2 I A!!/^ ^ lni^,[A,A2AMlD]+ E 



i=l 

ReS^2=^2(^1^2^3^4(^?)) 



poles{i} ^'^i 



(4.9) 



The first term encodes all the information from the boundary of the /i contour integral. 
The Inf operation therefore takes the form of a polynomial in /x which is cut off at some 
maximum power, 

ln^^.[f{^i)] = Y,cul?^ (4.10) 

fe=0 

The second term contains information about the pentagon coefficients since it has an extra 
propagator in //^ we can identify these terms with coefficients in the basis by performing a 
partial fractioning in //^ and comparing with equation ( 3.1C| ): 

Res^2=^2(AiA2A3^4(/l)) _ 
^2 - u2 



^2Res^2=^2(Ai^2^3^4(^r)) ReS^2=^2(AiA2^3^4(/"l)) 



(4.11) 



The first term in the last equation gives the contribution to the D + 2 dimensional pentagon 
integral whereas the second term combines with the ijP component of the boundary term 
to form the coefficient of the = 4 — 2e dimensional box integral. It is true in any case 
that neither of these terms contribute to the rational part and therefore we can extract the 
information from standard four dimensional cuts. Explicitly we find, 

^_ >.-R.V.;(A.A.A3A.)._„^_ (4.12) 
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which matches up with the vanishing of the D + 2 dimensional pentagons from equation 

dnoD. 

The rational contribution is therefore found by looking at the behaviour of the product 
of the four tree amplitudes at large values of /U, 

4 

o-=± A:=0 

^ = ^ E Inf^2[^i^2^3^4(rr)]U^, (4.14) 



<T=± 



where we have remembered that the quadruple cut of the scalar box integral is — i(47r)^/^. 
In the final formula, the Inf^2 operation has been restricted to the coefficient of the ii^ 
term of the polynomial. We discuss the maximum possible power of fi^ appearing in eq. 
( 4.13 ) in Appendix It is straightforward to look at the limit ^ — > oo by performing a 
Taylor expansion to give analytic expressions. 



4.2 Triangle Coefficients 




Figure 2: Momentum conventions for the triple cut; all momenta are outgoing and the loop 
momentum flows clockwise. 

We next consider the triple cut integrals^ which, with the new dependence on /x, can 
be written as, 



d%\{5{lj - ^?)A^A2A^ 

i=l 

I I dtJ,(lnf^.[lnU[AiA2As]]+Y,lnf^^ [Res,=, (^lA^Ag; 



t-ti 



^ Res„2=„2 Inft AiA2^3 ^ Res^2=^2(Rest=t,(AiA2^3)) 



(4.15) 



^We have suppressed the factors of i and tt^ in the triple and double cuts, these are implicitly put back 
in when writing the formulae for the coefficients. 
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Defining Kl 3 analogously to equation (|4.3| ), we choose a loop momentum basis, 



h = aKl + bK\ + t\Kl)[K\\ + 
which ensures that there all integrals over and l/t" vanish [28], 



/ 



dtjfe = 0, 



dtjt— = 0. 



(4.16) 



(4.17) 



The second and fourth terms can be directly associated with the previously calculated 
scalar box coefficients. This is due the presence of an additional propagator term and the 
lack of t-dependence in the numerator. We can also show that the third term does not 
contribute to rational part by applying partial fractioning in an analogous way to the box 
case. Of the two terms obtained, one will be associated with the 6 — 2e dimensional boxes 
and the other with the 4 — 2e triangle coefficient which is cut-constructible. 

Therefore the leading /x^ dependence of the triangle functions can be completely de- 
termined through the boundary behaviour: 



(4.18) 



One must also sum over the two solutions, a, for the loop momentum. For the massless 
case when 5i 7^ and ^ Qit is sufficient to sum over the two solutions for 713. However, 
in general it is necessary to sum over the solution li given above in eq. (P^ ) and conjugate 
solutions given by. 



II = aKl + hK\ + t\K\)[Kl\ + 



13 



713* 



\Kl)[K 



Ih 



(4.19) 



evaluated at a fixed value of 713. This applies equally well to the one and two mass triangles 
when one or both of Si, 53 vanishes and there is only a single solution for 713. As for the 
box case, we postpone justification of the form of eq. ( 4.18| ) to Appendix 



4.3 Bubble Coefficients 





K3 



Figure 3: Pure Bubble and Triangle terms contributing to the bubble coefficients. 
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Finally, we find the double cut integral using the following basis for the loop momen- 



tum: 



where 



h = yKl + 



Siil-y) 
7 



X + t\K\)[x\ + {y{l-y)Si-i?) 



\X)[K\\ 
It 



7 

The cut integral can then be decomposed into, 

2 

•' i=l 

= J d-^'fi J dtdy Jt,y Inf^2[Inft[Infj,[AiA2]]] 
+ ^Inf,. Lfj5ta(^ll+^lnf^, 

, Y^T f rRest=t (ReSj/=y,(Ai^2))l , v;^ 

+ > Inf,,2 — r + > 

^ ^ {y-yi){t-tj) ^ 



7 = 2(i^i-x). 



(4.20) 



(4.21) 



Rest=t,(Inf^[AiA2]) 
t-ti 

ReV=/.f(Inft[Infy[Ai^2]]) 
/^^ - A 



Res^2=^2(Inft[ReSy=2/,(^iA2)]) Res^2=^2(Rest=t,(Infy[AiA2])) 



{iJ^ - tJ^){y - yi) 

Res^2=^2 (Rest=t^. (Resy=y, [^1^2])) 
{f^'^-f^l)it-tj)iy-yi) 



(Ai2-/.2)(i_i.) 



(4.22) 



The last four terms contain residues in // and as such cannot contribute to the rational 
terms. Of these terms, the pure fi'^ — ^if pole will also have a contribution to the 11 = 
(or cut-constructible) bubble coefficient. The fourth term, in which the numerator is 
completely independent of y and t, can also be discarded as it has two propagator terms 
and can therefore only come from scalar box integrals. The remaining terms do have 
contributions to the bubble coefficient, and as described by Forde [28], can be determined 
by computing the Inf expansions for both y and t for both the double cut and triple cuts 
with non- vanishing integrals over t. The additional ji dependence is again determined by 
looking at the large ^ behaviour: 



^bub[2] ^ _iinf^,inf4lnf^[Ai^2(ri(2/,t,M')]]]U^to,j/^ 
<^tri(i.3)[2] = -1 ^ Inf^2lnf,[Al^2^3'''(ri(ya,i,/x')]]U^t»^T,. 



(4.23) 
(4.24) 



(T=± 



with the full bubble coefficient being a sum of the pure bubble and triangle subtraction 
terms: 



^ (7bub,[2] ^ ^ ^tri(fir3)[2] 
{^3} 



(4.25) 
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The functions Tj and Yi have been computed recently in references [28, 30] for arbitrary 
kinematics. Explicitly with a uniform internal mass we have, 



y„ = l y,= l Vi = l(l-|). (4.26) 

Solving the additional on-shell constraint, (Zi + K^Y — triangle subtraction 
terms gives the two solutions for y as, 

^^^(W^i^, (4.27) 

where 

C2 = S^{X\K^\K\1 (4.28) 

Ci = ^t{K\\K^\K\] - S^t{x\K^\x] + S^{x\K^\K\l (4.29) 

Co = -it^{K\\K^\x] - ^?{x\Kz\K\]+^tS^ + tS^{x\Kz\x]■ (4-30) 

The non- vanishing integrals over t are given by: 

_ s,{xmK\] 

2^A ' ^^-^^^ 

n = -'-^^^0^ (SiSs + K, . KsS,) , (4.32) 

n = - ^^j^af r (l5>gi>g| + 30Ki • ir3.5?.53 + ll(ifi • K^fSf 



+ SSfSs + IGn^SfAj , (4.33) 

where A = (i^i • K^ f - S1S3. 

5. Rational contributions to gluon amplitudes 

Integral coefficients in which a single particle type circulates in the loop are particularly 
well suited to the method described in the previous section. Such calculations apply to 
give the rational terms of the all-gluon amplitudes. The supersymmetric decomposition of 
such amplitudes is given by, 

= A^=^ - 4A^=' + + AT; [a^=' - 4^]) . (5.1) 

Since supersymmetric amplitudes are cut constructible in four dimensions the rational 
terms in such an amplitude only appear in An': 

41 = CW+K. (5.2) 
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J + 1 



J + 1 i 



' 

; ; fx XX 

; + 1 I + 1 

Figure 4: Cut diagrams for the integral coefficients contributing to the rational parts of the n-gluon 
amplitudes. 

We can then compute the rational parts of any scalar amplitude by introducing an effective 
mass and evaluating the integral coefficients in four-dimensions with the tree amplitudes 
for massive scalars [19, 69, 70]. 

The rational contribution, written in terms of the massive scalar basis, is therefore, 

^ n i—3 i—2 i—1 

i=l j=i+l k=j+l l=k+l 
^ n i—2 i—1 

i=l j=i+l k=j+l 
i=l j=i+2 

as depicted in figure ^. The arguments of the integral coefficients are considered to be taken 
mod(n), which accounts for the extra factors in front of the coefficients due to an over- 
counting in the summations. Our notation also suppresses the momentum flowing from 
one of the vertices of the integral since it can be inferred from momentum conservation. 

5.1 Analytic expressions for the four-point amplitude 

As an explicit example we present an analytic computation of the four gluon amplitude 
which has been considered previously with similar unitarity constructions [18, 45]. Here we 
only require the three- and four-point tree level amplitudes with a pair of massive scalars 
[69, 70]: 

A,{1,, 2-, 3.) = ifflM, ^3(ls, 2+ 3,) = .ffl^, (5.4) 

^.(1..2^3^4,)=i^^. A.(l,.2^3-,4,) = i||||J. ,5.5, 

where /i is the mass of the scalar particles. We then construct the integrand for the 
quadruple cut using K4 = p4,Ki = pi, K2 = P2 which leads to a simplified solution of the 
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general on-shell constraints from eq. 

a = 0, 6 = 0, d = — —. (5.6) 

CS41 

The fact that we only require information about the leading behaviour in ^u^ means it is 
sufficient to look at the leading behaviour of the coefficient c: 



In fact, since we have a product of four tree amplitudes, we find the coefficient is only a 
function of and so there are no square roots appearing in the final values: 

cW(l+,2+,3+,4+) = »l, (5.8) 

4 V , > , 7 (12)(34) ^ ' 

^4 (1 '2 ,3 ,4 )-- (24) (34) ' ^^'^^ 

cM(l-,2-,3+4+) = ^^MM, (5.10) 
4 V > > > y (34)[21] ^ ^ 

cWfl- 2+ 3- 2^(12) (34) [42]^ 

^4 (1 ,2 ,3 ,4 )- ^24)2 [21] [43] ' ^^"^^^ 

The triangle coefficients are particularly straightforward to evaluate since the Inf^ operation 
is equivalent to performing a Taylor expansion around t = oo. It is helpful to choose the 
two massless legs of each triple cut to form the basis since this simplifies the form of the 
analytic expression. Let us be slightly more explicit by giving the specific details for the 
computation of the C3;i2(l~, 2+, 3"'', 4"*") coefficient. In this channel it is easiest to choose 
Ki = p3, K2 = Pa since both Si and S2 will vanish and the two solutions, one each from 
eq. (|1^) and ( p^ , are: 



/i=t|3)[4|-i^|4)[3|, (5.12) 
/i* = t|4)[3|--^|3)[4|. (5.13) 
The product of tree amplitudes is trivial to write down: 

4+ , U r / 1- 2+ / U f / 3+ / 1 - g. (l|^2|2]^(3|/i|4](4|/i|3] 
^i(-/i,4 ,/2)^2(-/2,l ,2 ,/3)^3(-Z3,3 M)-2i (i2)(34)2[12](l|/2|2] • ^^-^^^ 

We are then left to insert the two solutions for li and sum over the two values according 



to equation ( [5.751) . 



Int,[Inf,[^,^2-43(;i)ll|„.,. = ■ (5.15) 
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So, after applying momentum conservation, we find the final result: 



C7gl2(l-,2+,3+,4+) 



I s 



12 



2si4) [32] 



(23) (24) (34) [21] [31]' ^^'"^^^ 

Following the same procedure for the other channels and helicity configurations quickly 
yields: 

i(12)[32]2 



C7£^3(l-,2+,3+,4+) = -. 

cgjl-,2+,3+,4+) 



•'3:41 1 



(24)2 [21] ' 
^gl2[32][43] 
(23)(34)[31]2' 



,3+, 4+) 



i {4?, 



2slA (12) 



and for the — I h amplitude: 

CP}2(l-,2+,3-,4+) 

(7^3(1-, 2+,3-, 4+) 

c||4(l-,2+,3-,4+) 

cgli(l-,2+,3-,4+) 
with all other coefficients vanishing. 



(23) (24)2 (34) [41] 

2iSi2S24 



(24)2 [31]2' 
2zg23(13)(23) 
(24)3 [41] ' 

cgUl-,2+,3-,4+), 



a 



[2] 

3;23 



;i-,2+,3-,4^ 



(5.18) 
(5.19) 
(5.20) 

(5.21) 

(5.22) 

(5.23) 
(5.24) 



For the bubble coefficients, we minimise the number of triangle subtractions by choos- 
ing X = pi for Ki = pi + P2 and X = P2 for Ki = P2 + Pz- Again, this has the benefit 
of giving us results directly in terms of the external momenta and yielding simple analytic 
forms. The the individual solutions for y+ and y_ in the triangle subtractions can lead to 
some spurious denominators and square roots. It is, therefore, beneficial to perform the 
sum over the two solutions algebraically. By considering the boundary behaviour at large 
t we can write the solution for y± as: 



y± = ai,±i + Oi2,± + -a3,± + ( ^ 



oii,± 



a2,± 



a3,± 



^{Kl\K3\Kl]-Si{x\Ks\x]±V^ 
2Si{x\K3\K\] 
a^{Si + 52) 




a 



Q:(x|A';,|A7](5i - l/rO v^v5i(x|A'3|A7](5i + S2) 



a 



a = {j{kI\K3\kI] + Si{x\Ks\x]f - 4f 



(5.25) 
(5.26) 
(5.27) 

(5.28) 
(5.29) 



The coefficient is then left as a function of Oj ± and by using sets of identities such as. 



1 



± 



± 



(7(i^^|i^3|i^i^]-gi(x|i^3|x])' + « 
4Sl{x\K,\K[]^ 



(5.30) 
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we are able to identify terms free from square roots after the summation is performed. 
Such a procedure is systematic and has been performed using symbohc manipulation in 
FORM [71]. Therefore, after a little algebra, we quickly arrive at the expressions: 







q+ /| + ^ _ 

1 "J ) 4 ) — 


2i(si3-S23)(13)(14) 
(12) (23) (24) (34)2 [21] 




^2;23(1 


',2+ 


,3+, 4+) = 


2iS24 (S24-S12)(12)[32]2 

.33(24)2[21] 


(5.32) 


C2;23(l 


,2- 


,3+, 4+) = 


2i(2si2-3s23)(12)2[41] 
3(14) (23)3 [21]2 [32] 


(5.33) 


'-"2;12'.-'- 


-,2+ 


,3-, 4+) = 


2i(5si2 + 2s23)(13)' 
3.2,(24)2 


(5.34) 


^2;23(1 


-,2+ 


,3-, 4+) = 


2i(2si2 + 5s23)(13)2 
3.23(24)2 


(5.35) 



with all other bubble coefficients evaluating to zero. 

The resulting expressions for i?| agree numerically with the known analytic results 
[45]. The form of the rational term is slightly different, and in some cases less compact, 
than those obtained through on-shell recursion but include both direct recursive and cut 
completion terms. Here we have the advantage that we do not need to calculate any 
completion terms for the cut-constructible parts and there are no problems associated with 
the factorisation in complex momenta. 



5.2 The five- and six-point all-plus amplitudes 

In this section we demonstrate how the technique applies to higher point amplitudes in a 
straightforward way. The finite helicity configuration with all gluons carrying positive helic- 
ity is a particularly simple example since there are no triangle or bubble type contributions. 
We only need a single five-point tree amplitude to evaluate these all-plus configurations up 
to the six-point level: 



,(0)n 0+ 0+ . N_ »/^^[2|l(2 + 3)|4] 

(23)(34)(2|1|2](4|5|4] 



4^\i.,2+,3+,4+,5.) = (5.36) 



We can again make use of eq. (|5.7|) to quickly determine the leading dependence of 
the quadruple cut which results in: 

i?f(l+,2+,3+,4+,5+) = -^cW2(l+,2+,3+,4+,5+) + cyclic perms. (5.37) 

ci'L(l^2^3^4^5-) = «!HM (5.38) 

^'^2^ , , , , ; (12)tr5(4,l,5,3) ^ ' 

where 

tr5(l,2,3,4) = (1|234|1] - (1|432|1] (5.39) 
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For the six point amplitude we have three independent coefficients which also follow from 
a similar procedure: 



i?f(l+,2+,3+,4+,5+,6+) = 

- ^C^23(1^2+,3^4^5^6+) - icW,|3,(l+,2+,3+,4+,5+,6+) 

- ^aS:;l2|45(l^' 2+,3+,4+,5+, 6+) + cyclic perms. (5.40) 

^[4] „+ .+ _ 2z (.45(6|1 + 2|3][51][64]^ - 546(5|1 + 2|3][54]^[61]) [56] 

C,.^,,,[l ,2 ,3 ,4 ,5 ,6 )- (12)(23)tr5(5,4,6,l)tr5(5,4,6,3) 

(5.41) 

Cfl 3,(1^2- 3-4- 5- 6+) = 2^(511 + 2|6](6|1 + 2|5][12][43][65F 
4;i2|34^ ' ' ' ' ' ^ (12)(34)tr5(5,2,6,l)tr5(5,4,6,3) ^ ' 

[4] + + + + + + _ 2z((3|l + 2|3](6|l + 2|6] 36^12) [12][54][63p 

^4;i2|45(l '2 ,3 ,4 ,5 ,6 )- (I2)(45)tr5(2,3,6,l)tr5(5,3,6,4) ^^"^^^ 

These expressions all agree with the previous analytic expressions [35, 45]. 



5.3 The five-point MHV amplitude 

As a more involved example we present analytic expressions for the five-point MHV config- 
uration i^Kl"*", 2+, 3-, 4~, 5~). This amplitude has been derived previously using a string 
based analysis [72] and more recently using on-shell recursion relations [37] 

i?f (1^ 2^ 3^ 4^ 5-) = ( - (,^,)(1+, 2^ 3+ 4^ 5-) - icg^^^^^ 

fe=i ^ 

_1^[2] /-,+ 9+ o+ A- ^k{k+l) „[2] „+ „+ ^.-A 

2^3;it(ik+l)|(fe+2)(it+3)V-^ '^^ ) g '-"2;it(jfc+l)'^-^ )} 

(5.44) 

= -^cw,.^-5cgou>-^:^cg,*«,(l^2^3^4-,5-) (5.45) 



g 4;total 2 3;total / y g 
fe=l 

The computation of box contributions follows in the same way as for the all-plus configu- 
ration though the lack of symmetry means each coefficient must be computed separately. 
The additional tree level amplitude required is given by [69, 70],^ 

4^1.,2-,3^4^5.)- ^(2|1(3 + 4)5|4]^ .,^[34]3 



(23) (34) (2| 1 12] (4|5|4] [2|(3 + 4)5|4] S234 [23] [2|(3 + 4)5|4] ' 

(5.46) 



*We note that one must compensate for the different norniahsation of the spinor products used in [69]. 



-17- 



The sum of all five box contributions then yields the following simple result: 



W;total-Z.W;fc(fc+l)U ,^ ,^ ,4 ,b J " ^^^^^ 3, 3, 4) V 

(14)2(25)^[53]^[21]3 _ 2(14) (25) (45) [31] [53] [21]^ (45)^ [31]^ [52] [21] (45) [31] [32] [21] 
(12)2 (34) [43] [51] [52] (12) (34) [43] [51] ^ (34) [43] [51] ^ 

(45)2 [43] [53] [21] 2(24) (35) (45) [31] [32]^ [41] (45)2[31]2[32][42] 
+ (15)(23)[43][51] ^ (15) [43] [51] 

(24)2(35)2[32]3[41]2 (45)2[32][41][51] \ 

(15) (23)2 [42] [43] [51] (23)[54] J' ^ 



while for the sum of the ten triangles we have, 



!] _ 2i(24)[32]2[21]3 i(23)[32]3[21]3 i[43][21]3 

t^t'^i ~ ~(2|5 + 1|2]2[42][51][52] " (2|5 + 1|2]2[42]2[51][52] " (23)[42]2[51][54]i 

i(12)[32]3[21]3 2i(24)2(25)[32][21]2 2i(25)[32]3[21]2 

(2|5 + 1|2]2[42][43][52]2 " (12)(23)(2|5 + 1|2]2[51] " (2|5 + 1|2]2[42][43][52] 

z(15)[53]=^[21]2 2i[31][21]2 i(24)2(25)2[32][52][21] 

(12)2[43][51][52]2[54] ^ (23) [42] [51] [54] (12)2(23)(2|5 + 1|2]2[51] 

2i(24) (25)2 [32]2 [21] i(24)2 (25)2 [32] [42] [21] i(45)[53]2[21] 
(12)(23)(2|5 + 1|2]2[43] ~ (12)(23)2(2|5 + 1|2]2[43] ^ (12)2 [51] [52] [54] 
i(15)[32][53]=^[21] i[31]^[21] i(45) [32] [41]2 i[31]^[32] 

(12)2[43][52]2[54] ^ (23) [43] [51] [54] ^ (23)2 [42] [43] [54] ^ (12) [43] [51] [54] 

i(34)[31][32][41]2 2^[31][32]2 i[32]^[5l] 

(23)2 [42] [43] [51] [54] ^ (12) [43] [52] [54] ( 12) [43] [52] 2 [54] ' ^ 



The bubble coefficients are more complicated for higher point amplitudes due to the ap- 
pearance of additional poles in the triangle subtraction terms. In order to give compact 
expressions it is convenient to leave the sum over y± unexpanded. We choose a value of 
X = Pk for the Ki = pk +pk+i channel which ensures there are a maximum of two triangle 
subtraction terms. Furthermore it is useful to note that to evaluate the coefficient at the 
boundary of the /x contour it is sufficient to use the following values for the non-vanishing 
integrals: 



Yo = l Y,=0 ^3 = -^ (5.49) 
= 2^A = ^ = 3^3A2 (^-^0) 



The justification for this can be seen be examining the expanded forms of the Inf operations 
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given in appendix 



Cn2(l+,2+,3+,4-,5-) = 2^([32][51] -3[31][52])[32] 



(12)2[21][43][52]2[54] 

(15)((24) -a};f(14))2(a};^^(15) - (25))[21]([31] + a^Jf [32])2[52] 

— 2i y ■ ■ , (5.51) 

.ti D33(aSM,l + 2)(12)2(34)(5|l + 2|5]2[43][51] 

^[2] ,^ _ 2i([31][42]-3[21][43])[31]^ 

02.23ii ,^ ,4 ,5 j- ^23)2[32][43]2[51][54] 

^ («^:^"(24) - (34))((35) - a^:yM25))^(2|5 + 1|3]([21] + ag^M31])2[32] 
^ti Z?55(«?;r>2,2 + 3)(15)(23)2(4|5 + l|4]2[42][51] ' ^ 



n+ 3+ 4- = 2i(24)(25)^(35) 10^(25)^(45) _ 4i(24)[21][41]2 

2;34^ ' '"^ ' ^ 3(12) (15) (23)4 [43] 3(12) (15) (23)3 [43] (23)3 [42] [43] [51] [54] 

4i[21][31][41] 2i[21][32][41]2 2i[21] [41]^ [53] 



(23)2 [42] [43] [51] [54] (23)2[42]2[43][51][54] (23)2[42][43][51][54]2 

■(ag-(35) - (45))(3|1 + 2|4][21][43](a?;r)^ a?;^^I)5i(«?f , 3, 3 + 4)2(3|4|2][42](ag^ )3 



cr=± 



Z)ii(afr,3,3 + 4)(12)(5|1 + 2|5]2[53] ' I)55(«?:f , 3, 3 + 4)(15)(2|3 + 4|2]2[51] 



^ 8^51 ,3,3 + 4)2 (23)2 (34) [42]^ ( [32] + ag^ [42] ) {a^f f 



+ 



3^)55(a?;^^3, 3 + 4)s34(15)(2|3 + 4|2]4[51] 
2D51 jatf , 3, 3 + 4)1)^1 (g^;^^ , g^;^^ , 3, 3 + 4) (34) [42] ( [32] + a'.'^J [42] ) (g^'g^ )3 

1)55(4?, 3,3 + 4)(15)(2|3 + 4|2]2[51] 
D5i(«S'^3,3 + 4)2Z)^5(g^;?,g^:?,3,3 + 4)(34)[42]([32]+g^:?[42])(g^:?)3 

I)55(a?;?,3,3 + 4)2(15)(2|3 + 4|2]2[51] 
gg^D5i(«^r,3,3 + 4)2.34[42]([32]+g;;?[42])(ga^)2 

Z)55(«?;?,3,3 + 4)(15)(2|3 + 4|2]2[43][51] 
2D5i(a?;r, 3, 3 + 4)2 (2«J?«?;^^.34 + l) [42] ([32] + [42])g?;? 

I)55(a?;?,3,3 + 4)(15)(2|3 + 4|2]2[43][51] " 



c|l5(l+,2+,3+,4-,5-) = 0, (5.54) 
cgi(l+,2+,3+,4-,5-) = cg4(3+,2+,l+,5-,4-), (5.55) 
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where we define the following functions, 

k;Ks _ (fc|fc + l\k]{k + l\K3\k + 1] - Skik+i){k\Ks\k] T a 

"^•^ " 2s,^,^,){k\Ks\k + l] ' ^'-^'^ 

«3±'=T— (fc|^3|A; + l], (5.57) 
' a 

D^y{a,p,q) = -a^{xp)[q^y] +a(^{x\q^\y] - i^\p\y]j + {x(i')[py]^ (5-58) 

D',y{a,p,p,q) = - + 2a(3^ {xp)[q^y] + /? (^{x\q^\y] - ^(^IpIv]) > (5-59) 



-,2 



where q = q p. 

2p-q 

The function a is given in equation ( |5.29| ) . These expressions have been checked numerically 
against the known results of [72]. 



5.4 The six-point — I i h amplitude 



As a final analytic example we consider the NMHV six-point amplitude with alternating 
helicities, i?g(l~, 2+, 3~, 4+, 5~, 6"^). This amplitude has been considered in previous an- 
alytic studies based on Feynman diagrams [73] as well as more recent numerical unitarity 
approaches [31, 53]. Here the analytic expression can be reduced to a set of seven inde- 
pendent coefficients through the symmetries of the external helicities. We first define three 
transformation identities: 

ai : (1, 2,3,4,5,6) ^ {1 + 1,1 + 2,1 + 3,1 + 4,1 + 5,1 + 6) (5.60) 
a] : (1,2,3,4,5,6) ^ {I + 1,1 + 2,1 + 3,1 + 4,1 + 5,1 + 6)\q^[^ (5.61) 
P : (1, 2, 3, 4, 5, 6) ^ (6, 5, 4, 3, 2, 1) (5.62) 

The rational contribution to this amplitude can then be written as, 
i?f(l-,2+,3-,4+,5-,6+) = -i Yl U (^823) + «fc (c^SLlsJ 

fc=0,2,4 \ 

+ Ok (C'il2|45) + "I+l (^4;123) + (^4;12|34) + "fc+l (^ILks) 

3ak (^C'3;1234) + 3afc (C'3;1232;45) + ^^^fc+l (^3;123a;45 ; 123x45)) 

+ '^"I (/^ (^3;123x45)) + (C'3;12z34x56) + (C'3;l2x34x56) 

+ Sk+l,k+2ak (C'2^12) + ■5fc+2,A:+3afc+5 (C'2^12)) + Sk+l±+2,k+3C'k 

The procedure for computing analytic forms for these coefficients has been automated with 
the help of symbolic manipulations in FORM [71]. In particular the procedure does not 
require any prior algebraic manipulation of the tree level amplitude although the final 



-20- 



form of the coefficient will depend on the form of the integrand. We have found it use- 
ful to choose a basis for the loop momenta that reflects the symmetry of the coefficient. 
Therefore adjacent massless legs were chosen where possible. However for the 2-mass easy 
box configuration it is much simpler to choose the two, non-adjacent, massless legs. The 
solution to the on-shell constraints in this case can be easily obtained from ( ^.51 ) using an 
exchange of K2 ^ K^. 



cg23(l-,2+,3-,4+,5-,6+) 



tr5(r/i,2, 5, 4, 6)2tr5(7?3,2, 5, 4, 6)^ [64] (13)4(5|46|5)[64] \ 



(5.64) 



(5.65) 



(4|6|5](6|4|5][21][32]tr5(5, 1, 6, 4)tr5(5, 3, 6, 4) si23(12)(23)(46) J 
[4] . + + ^ 2i tr5(r?i,2, 5,1 + 2, 6)^tr5(r?3,4, 5,1 + 2, 6)^ 

cS2|45(l"'2+,3-,4+,5-,6+) = cg2|34(l-,2+,5-,4+,3-,6+) (5.66) 

where we have defined complex momenta -q^- = and ffi'j^^k ~ ^(^l^^il^)- 

Following a similar automated procedure expressions for the triangle and bubble con- 
tributions have also been generated and have been checked numerically against the known 
results. However, since the expressions for these remaining coefficients are quite lengthy 
and not particularly illuminating, we include a Mathematica input file, R6mpmpmpCoef f s .m, 
containing the coefficients which follows the notation of the SSM package [63]. To present 
more compact forms ofthe triangle subtractions terms it was necessary to use a slightly 
expanded basis functions than was used in the five-point example: 

D^,yAa,p,q) = -a\xp)[q'\y\z) + a (^{x\q'y\z) - ^^|^^) + {^Q%\y\z), (5-67) 
D'ia,p,q) = -a^[xq']{p\y\z] +a( [x\q'y\z] - ^^^^ + [xp]{q'\y\z], (5.68) 



<.(a,/3,P,9) = - + {xp)[q'\y\z)+(3(^{x\q'y\z) - fi^M^ ) , (5.69) 

<.(«,/3,P,9) = - + 20/3) {xq'){p\y\z]+p(^[x\q'y\z] - ) , (5.70) 



with all other definitions given in equation ( ^.601) . The bubble coefficients are then given 
by: 

/^[2] _ .^[2],bub 1 ^[2],-fC2=P6 , ^[2],-fs:2=P5+P6 , ^[2],E-2=P4+P5+P6 (r,7^\ 

'-'2;12 — ~**-'2;12 ~ 9 '-'2;12 + '-"2;12 + '-^2;12 l^-'^^j 



*-^2;123l-'^lJ — ~^*-"2;123 



- I E CS!{2T'%KI) + Cat=^^+^«(i^J) + cf^,f='^'-'^-Kl) (5.72) 
where X = Pi has been chosen for the basis in both cases. 
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Since the procedure for analytic extraction has been automated it has also been possible 
to generate unsimplified expressions for all other configurations with up to six gluons which 
have been checked numerically against existing results. These expressions are available from 
the author on request. 



5.5 Remaining helicity configurations 



Since the extraction is purely algebraic is also well suited numerical extraction. One possible 

4 



way to do this is to use the discrete Fourier projections to find the value of c|^^ . From eq 



(4.14) the box contribution can then be written as, 



= ^ E E ^^1^2^3^4(/"?,//fc), (5.73) 
^Pl" a=± k=0 

^J'k = /^ooexp (27ri^^ , 

where ^oo is some large constant and + 1 is the number of points around the circle 
of integration, on which the residue is evaluated. It is straightforward to write down 
equivalent expressions for the triangle, 

^ C'P = -^E E E -,A^A,A,(ll,t,^^?,\ (5.74) 

^J'k = fJ-oo exp (27ri^^ , U = t^o exp (^nij^^ , (5.75) 

and bubble terms, 

PM-lpt-lPy-l 2 

C^2"'''^ = — E E EE^^Mh{y,,t,-i,4)) (5.76) 

= - ^7 E E E E -^MAMWiy^^t,,.,l)) (5.77) 



/ioo exp (27ri^^ ti-^ = C exp {^T^iy^ Vq = Voo exp {^m-^ 



Since the locations of the boundaries depend upon the order in which we take t and y large, 
we must take a = —1 for the bubble terms and a = 1 in the triangle subtraction terms. 

The method of Fourier projections has already been shown to be extremely efficient 
in both the BlackHat code [31] and with the OPP technique [25]. In general one can 
arbitrarily increase the numerical accuracy by increasing the size of the circle, fioo, and the 
number of points > 4. We note that, in order to find a stable numerical implementation 
at fixed precision, it may well be necessary to make explicit subtractions of the poles from 
pentagon integrals [21, 26, 53]. This could be done in an analogous way to the treatment of 
the cut-constructible contributions in the BlackHat code by computing values of /x^ which 
put the additional propagator on-shell [31]. 
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We have tested our method for the direct extraction of the rational terms for the 
remaining one-loop gluon helicity amplitudes with up to six external legs using a numerical 
approach with discrete Fourier transforms. In order to compare with the known results 
it was sufficient to proceed without the subtraction of higher order residues in place of 
considering a large radius for the integration contour. The required tree level amplitudes 
have been generated using on-shell recursion relations literature [69, 70]. The rational terms 
have been checked against the known analytic results [37, 38, 40] and the numerical results 
of references [52-54]. We have also checked these results against expressions generated 
analytically through an automated implementation of the direct extraction method. In 
order to achieve numerical precision of at least 10~^ it was necessary to use a radius for 
each contour integral of 10^ (i.e. the parameters fJ.oo,too,yoo) and in some cases evaluate 
at a large number of points around the circle (i.e. p^,Pt,Py ~ 10). This was achieved 
using symbolic manipulation in FORM [71] and numerical procedures in Maple. In order 
to obtain a more stable implementation, subtraction of the residues of high point functions 
from the complex plane [23, 31] would allow setting all the radii to 1 and minimise the 
number of points required around the circle. It is expected that these improvements would 
considerably improve the speed of the algorithm but we leave such an implementation for 
future study. 



6. Rational contributions to quark amplitudes 

As an example of how the method can apply equally well to amplitudes with external 
fermions we re-compute the rational parts of the four-point process with a pair of massless 
quarks, gg — > qq. The n-point colour ordered amplitude can be written as: 

The supersymmetric decomposition ensures that the only rational contributions come from 
the scalar term, An'^ (just as in the gluon case), and a\^^ where the fermion line follows the 
shortest path through the loop (see reference [65] for the definitions of the left and right 
moving primitive amplitudes). 

Figure ^ shows the rational part of a\^^ written in terms of the massive integral basis. 

In the following calculation we show that it is possible to reconstruct the rational 
contributions to the left moving primitive amplitude by looking only at the transverse 
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Figure 5: Integral basis for the rational part of the leading colour contribution to the gg 
amplitude, i?4^'. 



polarisation of the D-dimensional gluon. This can also be interpreted as introducing a 
massive scalar in the loop in an analogous way to the gluon examples. We then simply 
interpret the D-dimensional internal fermion as four-dimensional massive fermion ^. 

We now construct the tree-level structures using both massive scalars and massive 
fermions so that each cut has a uniform internal mass, 



Vsilq, 2s, 3q) = iu{pi)v{p3,n) 
V3(lQ,25,3g) = iu{pi,fi)v{p3) 

where we can use the following representations for the fermion wave-functions: 



{r] T \{<^ + ^J-) 



v±iq,r]) 



The relevant four-point amplitudes are given by: 

m(pi)t;(p4,At)(l|3|2] 



A4(1„2+,35,4q) 
A4(1„2",35,4q) 
^4(1q,25,3+,4,) 
A4(1q,25,3-,4,) 



(2|3|2](12) 
iu{pi)v{p4,fj,){2\3\l] 

(2|3|2][21] 
m{pi,lj)v{p4)m3] 

(3|2|3](34) 

m{pi,l2)v{pMm 



(3|2|3][43] 



^4(1+,2s,35,4-) = i(4|2|l] 



+ 



1 



SAi 2(1|2| 



(6.4) 
(6.5) 



(6.6) 

(6.7) 
(6.8) 
(6.9) 
(6.10) 
(6.11) 



^We do not present a formal proof of this fact since in general there is an ambiguity in the definition of 
the higher dimensional Clifford algebra. The definition is scheme dependent but while working within the 
FDH scheme we simply note that there is no problem in treating the D-dimensional fermion as a massive 
one for the example presented here. Indeed D-dimensional massive fermions have been treated successfully 
within the FDH scheme in a recent paper [55] 
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The leading mass components of the seven integral coefficients, shown in figure |5[ are 
easily computed with the techniques used for the gluon amplitudes in section |5|. Firstly we 
note that for the leading colour primitive amplitude all box contributions are zero since 
the maximum power of /x that can appear in the product of the four three-point vertices 
is three. The first non-trivial coefficients are therefore the triangle coefficients which we 
present here for the three independent helicity configurations. The -|- + +— configuration 
is given by, 



2 (23) [24]' 

^3;23U, '-^ '^5^- 2(12) [24]' 

^[2] 2+ 3+ 4-^ = ^^^^)[^^] + ^(14)^ [12] 

3;34^^9'^ ' ' «^ 2(23)(13) (12)(34)(13)' 

[2] „+ ^ ^ (^14(13) (24) + {su - 4^24) (12) (34)) [32] 

L'3;4iUg>^ ' ' 5^ 2(12)(13)(14)(23)[41] 

[2] n+ 0+ Q+ ^-^ - ^[31](2si3 + g23) 



X] '^3;fc(fc+l) (^9 , 2+ , 3+ , 4 
k=l 

the -|- H configuration by. 



(12)(23)[41] 



^[2] . + + Q- 4~^ _ ^(34)^ 

L'3;12Ug>^ ' ' ^ " 2(12) (24)' 

[2] 2+ 3- 4-1- !M(23)(34^ , ^(34)2[21] 

'-^3;23l-'-g '-^ J — o/lo\2/0/l\ ^ 



2(12)2(24) (23) (24) [32] 



^3;34(lg^' 2^) 3 , 4^ ) — C^Jgllj)"! 2^, 3 ,4^ 
fl+ 2+ 3- 1-)- ^<14)(23)[21] 

L'3;4lUg,^ ' ' ^ " 2(12)2 [31] ' 

V"r[2] n+ 9+ q+ 4-^ - ^^13(34) [21] 

Z.^3;fc{fc+1)^^5'^ '"^ '*'?^" 523(24) [43] ' 

and finally the H 1 — configuration, 

rm n+ 2- 3+ 4- ^_ ^(24) [41] [43] 

'-"3;12Uq ; ^ ) -J , % J 



2(23) [42] 2 ' 

[2] _ ^[31] (g24([31][42] + [21][43]) - g23[21][43]) 

2(23)[21][32][42]2 



f^lA (-1+ o- 0+ 4 



4 



^3;34(l<j^'2 , 3^,4^ ) — C3.|2(lg^i 2 ,3^,4^), 

[2] . + + . _ ^(14) (23) [31] 

t-3;4lUg,^ ,3 , 4^- j - 2(13)2[21] ' 
[2] n+ o+ o+ .-^_^(12)[31] 



3 



X^'^3;fc(fc+1)(^9'^^'^^'^9^ S23[24] 

There is only one non-zero bubble coefficient for all the four point fermion amplitudes 
considered here. This is the 23 channel for the + -|- H — helicity configuration for which we 
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choose a basis defined by Ki = P2+P3 with x = Pi K\ = p^). This choice ensures that 
there is only one triangle subtraction term and the two integrands are: 

C2;23(l,,2 ,3 '4,)- (23)(2|Zi|2] Wi^Wm) ^ ^ ^ 

^2;23 (1,,2 ,3 ,4,)- ^23)(2|/i|2] " ^^"^^^ 

Since Tq = 0, the triangle subtraction term will vanish after substitution and expansion 
of the loop momentum. The second term in the pure bubble is also suppressed by an 
additional power of the loop momentum and therefore vanishes after taking the large y 
limit. The single remaining term is then gives the final value for the coefficient: 

Cl|3(1^2^3^4.) = -?«i (6.29) 

The rational contributions to the left moving primitive amplitudes for the process qqgg are 
therefore given by: 

, 4 

i?P(l„2,3,4,-) = --^Cgi(,^,)(l„2,3,4,-) - ^0^3(1,, 2, 3, 4,-) (6.30) 

k=l 

Summing all of these components together can be quickly shown to match the results for 
four point amplitudes [65, 74]: 

[L]. + + + . _ ^(41) [13] / 2^12 \ 
it;4 (1,,2 ,3 ,4^-j- 2(12)(23) 3s23y' ' ^ ' 

i?P(l+ 2+,3-,4r) = l4°)(l+ 2+,3-,4r), (6.32) 
i?fJ(l+ 2-,3+,4r) = -1 (1 + J^) Af (1+ 2-,3+,4r). (6.33) 
The tree amplitudes above are the well known MHV amplitudes: 

A(°)fl+ 2+ 3- 4-) - ^^34)3(31) 

[1^,2 ,3 '4,-)- ^^2)(23)(34)(41)' ^^"^^^ 

a(°Vi+ 2- 3+ 4-) - '^^^^'^^'^ (6 35) 

(1,,2 ,3 .4,-)- ^^2)(23)(34)(41)- ^^"^^^ 



7. Conclusions 



We have presented a general method for extracting rational contributions to gauge-theory 
scattering amplitudes using D-dimensional unitarity techniques [26, 53]. Exchanging the 
D-dimensional cuts for four dimensional massive cuts, we then used simple complex anal- 
ysis to express the D-dimcnsional integral coefficients as contour integrals over a complex 
mass, evaluated on a circle at infinity. Here we have gone beyond previous approaches 
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[21, 26, 53] by making maximum use of the complex behaviour of the corresponding contour 
integrals. We have shown that it is not necessary to compute the pentagon contributions 
explicitly, thereby reducing the computation of rational terms to the evaluation of gen- 
eralised massive cuts of box, triangle and bubble functions in the large mass-limit. The 
formalism leads directly to compact analytic expressions for the rational terms. 

To demonstrate the method we have computed analytic expressions for the rational 
parts of one-loop gluon amplitudes using a massive scalar loop. This has been tested by 
re-computing gluon helicity amplitudes with up to six external legs. We have also shown 
that the method can also be implemented numerically using a discrete Fourier transform. 
The implementation was sufficient to evaluate all helicity configurations with up to six 
external gluons but we leave a detailed analysis of accuracy and speed to future work. 

We have also presented a simple example of how the method can be used to compute 
amplitudes with massless external fermions, considering the leading colour primitive am- 
plitude for the process qq gg. Here we show that the rational terms can be extracted 
using tree amplitudes where a massive fermion couples to a massive scalar inside the loop. 

Since the procedure is based on looking at four-dimensional cuts it is reasonable to 
expect that one can obtain a similar level of computational speed in comparison to calcu- 
lations of the cut-constructible terms [31, 47], though further analysis of the subtraction 
terms is expected to be necessary. It does have the benefit over on-shell recursive tech- 
niques that it avoids problems of unknown factorisation properties in the complex plane 
and therefore applies generally to any process, in particular to those with internal masses. 

Although the present paper concentrates on applications to amplitudes with massless 
external particles the main part of the procedure applies equally well to the case of massive 
ones. In this case the bases for the loop momentum in each of the cuts would have to be 
modified slightly as described in reference [30]. 
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A. fj,"^ Dependence of the Integral Coefficients 

In this appendix we use a simple power counting argument to demonstrate the maximum 
power of that can appear in the n-point gluon amplitudes fits the format described 
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in section |^. Although the final result is assured for general amplitudes in any gauge 
theory the following analysis sheds light on the analytic expansions of the Inf expansions, 
showing that only a few terms in the expansion contribute to the final rational terms. 
This information can then be used for efficient computations of the boundary behaviour. 
For the form to hold we must show that the leading dependence on fi appearing in any 
massive box cut is //^ and //^ in any triangle or bubble. This property is shown to hold 
in both the OPP/GKM formalism [26, 53] using a simple argument from the dimension 
of possible tensor structures. It is this argument that allows the parametrisation of the 
cut in the GKM approach and so we see that it is related to the decompositions of the 



Inf^ components of the cut given in equations ( 4. 13 ), ( 4.18 ) and the analogous relation for 
the bubble coefficients. In the analytic approach of Britto et al. the ;U-dependence can be 
determined through explicit analysis of the spinor integrals [21]. Here we use a slightly 
different approach to justify these statements, appealing to properties of the tree level 
amplitudes. 

We begin by considering a general quadruple cut given by the product of four tree 
amplitudes. Schematically the loop momentum can be written as, 

l[IJ, ) = ai + a2C -\ (A.l) 

c c 



the solution for c is given by equation ( |4.5| ) from which it is clear that c scales as y /x^ . It 
is then clear that the loop momentum given above scales as \//? in the large fi limit. 



The next step is to consider the dependence of a generic tree amplitude on the loop 
momentum. For clarity we restrict ourselves to considering a single massive scalar loop 
although it should also be possible to apply a similar argument to more general cases. The 
tree level amplitude can be characterised by considering a massive scalar flowing through 
a general Feynman graph as shown in figure The gluon-scalar-scalar vertex is simply 
proportional to Z • and each of the propagators will scale as, 



The maximum dependence on / will come from graphs with the maximum number of 
vertices. Such graphs will have n vertices but will also have n — 1 propagators, from which 
we can deduce: 

Af\l)^^^ XAf\l). (A.2) 

Taken together with the scaling of the loop momentum, this is sufficient to show that the 
product of four tree amplitudes scales at most as fi^ in the limit fi ^ oo. 

For the triangle coefficients we have to consider the expansion in t in order to observe 
the expected behaviour. Here the loop momentum is given as 

l{fi^,t) = ai + a2t+'j- + ^^- (A.3) 
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P2 P3 Pn-2 Pn^l 

•\.i:i....Li:./~ 



I I + Ki^n 

Figure 6: The momentum flow of a massive scalar through a generic tree level Feynman graph. 

We have that in any triple cut the maximum power of t that can appear will be three, which 
follows from the argument given above for the scaling of the tree level amplitudes. More 
specifically we have maximum difference of three powers of the loop momenta between the 
numerator and denominator and so we can write the boundary behaviour of a general triple 
cut as: 

InVlInt,[C3ll = (A.4) 



n n—2m 



Pn{x) = E E *V'"^^n-fe,m (A.5) 



m=0 k=—n 



In order to encode all the mass dependence in the coefficients a, b we have first re-written 
all explicit mass dependence in the tree amplitudes as fj,^ = li ■ li. After expansion, the 
component of C3 is, 



J f 1] 2/ ai o&2,i ao,o&3,i , ^ 
Inft[C3]|iO = n ^ -2 + 2 



^1,002,1 


Q3,i 1 


bo,o J 




^1,0^0,0 


"0,0 


bio 



\ "0,0 "0,0 "0,0 

Ql,0^2,0 I r) Q0,0^1,0^2,0 00,0^3,0 , ^3,0 ^l,0Q2,u , -i,u -^," i.u / A «^ 

"0,0 "0,0 "0,0 "0,0 Oo,0 

which shows that the triangle coefficients have the expected behaviour. 

Finally we use a similar argument to prove that the bubble coefficients also scale at 
most as /x^. The loop momentum in this case is 

2 2 
l{li^,t,y) = aiy + +03?/ + a^t + + (A.7) 

We then write down the boundary behaviour of the most general double cut integrand as 
a ratio of two polynomials, 

Inf^[Inf,[Inf,[[C2]]] = (A.8) 

n n—2m n—k—2m 

P'nix) = E E E /^'"*'y'^2n-;,fc+n,m (A.9) 

m=0 k=—n l=max{0,—k—m) 
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Again, using the argument on the scahng of the tree-level amplitudes, we see that the 
polynomial in the numerator has two powers more than the one in the numerator. The 
series expansion of C2 then gives us the following expression for the pure bubble cut: 



T r FT r 111 i^ ( '^2.2,0 00.0,0^2,2,0 «!. 1,0 ^1,1,0 , ^0,0,0 &1, 1,0^ 
lnU[lniy[C2\]W\y^Y, = Y2 \ i : o : n h 



^0,0,0 



^0,0,0 



^0,0,0 



h 3 
00,0,0 



_l_ Y ( — "0.0.0 ^2,2,0 _ Qi,i,o ^1,1,0 _|_ ao,o,o ^1,1,0' 



^0,0,0 ^0,0,0^ 



^0,0,0 



^0,0,0 



2 /^r, Q0,0,0 ^2,0,1 ^2,2,0 , Q0,0,0 fel, 1,0 ^3,1,1 02,0,1^2,2,0 , "4,2,1 00,0,0^4,2,1 Oi^i^Q ^3,1,1 



^o,o,o" 



^0,0,0 



^0,0,0^ 



^0,0,0 ^0,0,0^ 



^0,0,0^ 



03,1,1^1,1,0 02,2,0 ^2,0,1 , 02,0,1^1,1,0 ^ Oq^q^q ^j^i^Q ^2,0,1 , ^ Ol, 1,0 ^1,1,0 ^2,0,1 



1,0 ^0,0,0^ ^0,0,0^ ^0,0,0^ ^0,0,0^ 

00,0,0^4,2,0 Oi^i,o&3,l,0 Oi,o,0 ^3,2,0 02,0,0 ^2,2,0 02,1,0 ^2,1,0 02,2,0 ^2,0,0 



^0,0,0^ 



^0,0,0^ 



^0,0,0 



^0,0,0^ 



^0,0,0 



^0,0,0 



^0,0,0 



03,1,0 ^1,1,0 02,0,0 ^1,1,0^ 02,2,0 fei,o,o^ _ 03,2,0 fei,o,o 04,2,0 „ Oo,o,0 &1,0,0^&1,1,0^ 



^0,0,0^ 



^0,0,0 



^0,0,0 



&0,0,0 



&0,0,0 



&o,o,o" 



„ 00,0,0 ^1,0,0 ^3,2,0 r, 00,0,0 ^1,1,0 ^3,1,0 „ 01,1,0 ^1,0,0 ^2,1,0 „ Oi,o,0 ^1,0,0 ^2,2,0 



&0,0,0 



&0,0,0 



&0,0,0 



^0,0,0^ 



„ Ol,o,0 ^1,1,0 ^2,1,0 „ 00,0,0 ^2,2,0 ^2,0,0 „ Ol,l,o ^1,1,0 ^2,0,0 „ "0,0,0 ^1,0,0 ^2,2,0 



6n 



0,0,0 (^0,0,0 (^0,0,0 "0,0,0 

o Oo,0,0 ftl, 1,0^^2,0,0 ^00,0,0^1,0,0^1,1,0^2.1,0 o Oi l 61,0,0^^1.1,0 , r, 02,1,0 61,1,0 ^1, 0,0 

3 -_ 3 6 : 3 ■■ 3 : — ■ h 2 5 



60,0,0 (^0,0,0 
2 _ _ A. 2 



60,0,0"^ 



60,0,0"^ 



"0,0,0 



_l_ 00,0,0 62,1,0 _ 2 01,0,0 61,1,0 61,0,0 



60,0,0 



60,0,0 



(A.IO) 



Once we substitute the functions Yi from eq. ( f4.26| ) we see that the coefficient scales as 
For the triangle subtraction terms the parametrised loop momentum is the same as in the 
scalar triangle coefficient but this time we are interested in the coefficients of t^,t'^ and t: 



T f rr^Mii ( ^0,o\ , rr ( "l.O 090619 

Inf,[C, ]I*-T.=T3(— +T2 — 



^0,0 



^ f 2( 09,062,1 ^ 02,1 \ ^ 02,0 00,962,0 61,901,9 ^ 61909,0 

6g 60,9 / 69,0 6L 6^ 6^ , 



(A.ll) 



^,0 



^0,0 



^0,0 



which after we use equations ( [4.31 - 4.33 ) shows that the full bubble coefficient scales as 
expected. 
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